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Abstract 

We review the fundamental results of a new cosmological model, based on conformal gravity, 
and apply them to the analysis of the early data of the Pioneer anomaly. 

We show that our conformal cosmology can naturally explain the anomalous acceleration of the 
Pioneer 10 and 11 spacecraft, in terms of a local blueshift region extending around the solar system 
and therefore affecting the frequencies of the navigational radio signals exchanged between Earth 
and the spacecraft. On the contrary, conformal gravity corrections alone would not be able to 
account for dynamical effects of such magnitude to be capable of producing the observed Pioneer 
acceleration. 

By using our model, we explain the numerical coincidence between the value of the anomalous 
acceleration and the Hubble constant at the present epoch and also confirm our previous determi- 
nation of the cosmological parameters 7 ~ 10 -28 cm -1 and 5 ~ 10~ 4 — 10 -5 . New Pioneer data are 
expected to be publicly available in the near future, which might enable more precise evaluations 
of these parameters. 
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I. INTRODUCTION 

The Pioneer 10 and 11 spacecraft were launched in the early 1970s, to conduct explo- 
rations in the region of the solar system beyond the orbit of Mars and to perform close 
observations of Jupiter. They were also the first spacecraft to explore the outer solar system 
and to send back to Earth their navigational signals for almost thirty years (for a review see 
[l| and references therein). 

In recent years, the orbits of Pioneer 10 and 11 were reconstructed very accurately, by us- 
ing the original radio-metric Doppler tracking data, based on the signals exchanged between 
the spacecraft and NASA's terrestrial tracking stations. This reconstruction yielded a per- 
sistent discrepancy between the observed and predicted data, equivalent to an unexplained 
small acceleration of the spacecraft in the direction of the Sun. This effect is evidenced by 
measuring a small frequency shift (toward higher frequencies, i.e., a "blueshift") of the signal 
reaching us from the spacecraft. The nature of this anomalous acceleration or of the related 

n 

)lueshift remains unexplained; this effect has become known as the "Pioneer anomaly" ([2J, 

The importance of this effect, as well as of other known gravitational anomalies ([5], 
, is not related to how they affect the spacecraft navigation, since they all produce 
very small corrections to the orbits, but to the possibility that these anomalies might be an 
indication of new gravitational physics. In particular, several non-conventional explanations 
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of these effects have been proposed (see discussion in j^j], jj]) such as modifications of 
the law of gravity, or a modified inertia, as proposed by the Modified Newtonian Dynamics 
(MOND) theory, the existence of a dark matter halo around the Earth, or in the solar system, 
which might slightly alter the gravitational force acting on the spacecraft, and several others. 



In this line of reasoning, alternative gravitational t 
(CG), originally proposed by H. Weyl in 1918 ( 8|, 



reories such as Conformal Gravity 



10[) and revisited by P. Mannheim 



111 ]. 12], 13]), have provided a new framework for cosmological models, with the advantage 



of avoiding some of the most controversial elements of current standard cosmology, such as 
dark matter, dark energy, inflation, and others. 

Following the original CG, we have recently studied an alternative approach to these 
models which was named "kinematical conformal cosmology" Jjj], but that for brevity will 
be called Conformal Cosmology (CC) in the rest of this paper. This approach was based 
on the direct application of the conformal symmetry to the Universe, i.e., considering the 
possibility that a "stretching" of the spacetime fabric might be acting over cosmological 
scales. In a second part of this work 15] it was shown that this model can successfully fit 
type-la Supernovae data, without assuming the existence of dark matter or dark energy. 



In addition, a preliminary analysis also performed in our second paper [15] indicated 
that CC might be able to explain the existence of the Pioneer anomaly, since the observed 
blueshift of the spacecraft signal could be due to a region of cosmological blueshift surround- 
ing our solar system, which is naturally predicted by our model. A new comprehensive review 
of the Pioneer anomaly has recently been published [l|, together with more details of the 
^ioneer early data [16], thus prompting us to reconsider and improve our previous analysis 
15] , based on the conformal cosmology approach. 

In the next section we will briefly review our CC solutions, showing how a local blueshift 
region can naturally emerge, while in Sect. IIHI we will fit all current Pioneer data 1(| with 
our cosmological solutions and determine the values of the parameters in our model. Finally, 
in Sect. IIVI we will discuss our results and compare them to the existing physical limits of 
standard gravity in the solar system. 



3 



II. CONFORMAL COSMOLOGY 

In our first CC pa per [bjj ] we used as a starting point the line element originally derived by 
Mannheim-Kazanas as an exterior solution for a static, spherically symmetric source in 
conformal gravity theory, i.e., the analogue of the Schwarzschild exterior solution in general 
relativity: 

ds 2 = -B(r) c 2 dt 2 + J^y + r 2 d^ 2 (1) 
where dip 2 = dd 2 + sin 2 9 dcfi 2 in spherical coordinates and 

B( r ) = 1 - /3(2 ~ 3/37) - 3/3 7 + 7 r - kt 2 , (2) 

with the parameters (3 = (cm), 7 (cm -1 ), k (cm -2 ), where M is the mass of the (spher- 
ically symmetric) source and G is the gravitational constant. Conformal gravity introduces 
two new parameters 7 and k which are not present in standard general relativity, while the 
familiar Schwarzschild solution is recovered in the limit for 7, k — > 0, in the equations above. 

We then considered regions far away from matter distributions, thus ignoring the matter 
dependent (3 terms, and rewrote the last equation in a simplified form: 



B(r) = 1 + 7r — kt = 1 + + ^— + kj r = -g o(r), (3) 

2 

where the parameter k is linked to 7 and k, through k = — % — k and it is ultimately 
connected to the so-called trichotomy constant k (in bold) of a Robert son- Walker (RW) 
metric, defined as k = A = 0, ±1. This is also related to another fundamental aspect of CG: 
the existence of coordinate and conformal transformations connecting the static, spherically 
symmetric solution represented by Eqs. (OQ) and ([3]), with the classical Robertson- Walker 



metric (see details in }14|). 

It was precisely this connection between the two solutions which prompted us to consider 
the CG static, spherically symmetric solution as an alternative description of the standard 
cosmological evolution, based on the RW metric. In other words, the CG static solution 
might also contain information about the cosmological redshift, the expansion of the Uni- 
verse, etc., and constitute an alternative approach to cosmology. In particular, the CG 
expressions in Eqs. (j2J) or ([3]) contain a linear and a quadratic term, in the radial coordinate 
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r, which might yield considerable effects 1 at large distances, including a strong gravitational 
redshift which could be, at least in part, responsible for the observed cosmo logical redshift. 

Therefore, we postulated in Ref. [14| that the observed redshift is due to this gravitational 
effect, which influences the wavelength or frequency of a light signal emitted at time t and 
position r, and observed at the origin (r = 0) at the current time to, in the following way: 



_ R(0) _ X(r,t) _ u(0,t ) _ / -goo(0) _ 1 

R(r) A(0,t ) u(r,t) \j -g 00 (r) ^j l + ir+ ^ + k y, 

In the previous equation the redshift factor (1 + z) is related to the ratio of cosmic scale 
factors R, which simply depend on the radial distance r, in view of Eq. (J3]). Alternatively, 
to obtain a time dependent form of the cosmic scale factor, we considered that the radial 
distance r is associated with a look-back time (to — t), related to the time of travel of a light 
signal. Integrating the CG metric in Eqs. ([T]) and (j3J) along the null geodesic, we obtained 



14j: 



R(t ) 

1 + z = -— = (cos x ~ 5 sin x) ; k > (5) 
R(t) 

i + * = f| = (i-'x);* = o 

R(to) 

1 + z = = (cosh x ~ S sinh x) \ k < 0, 

for the three possible values of the parameter k. In the previous equation we preferred to 
use dimensionless quantities and parameters, defined as follows: 

, 'V 

X = V\k\c(t -t) ■ S = (6) 

2\/\k\ 



so that the fundamental parameters of our conformal cosmology are now expressed by 7 
(cm -1 ) and the dimensionless 5 (c is the speed of light in vacuum, assumed constant). 

In Fig. [T] we plot the results of Eq. (jSJ) in terms of the inverse ratio R(x)/R{0) = 
R(t)/R(t ) = 1/(1 + z) which describes better the cosmic evolution. The dimensionless 
quantity x = \/T^l c (^o ~ 0> 011 the horizontal axis, represents a look-back time, so that the 



1 By using the contribution of the linear term yr, the fiat galactic rotation curves were in fact explained by 
Mannheim ([l3|; 01) without the need of dark matter. 
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FIG. 1: R functions obtained from Eq. ([5]) are shown here for different values of k: k = —1 in 
red (solid), k = in green (dotted), and k = +1 in blue (dashed), and for a positive value of the 
parameter 5 ~ 0.6 (an unrealistically large value; our current value 5 = 6(to) will be shown to be 
positive and close to zero. 

universal evolution of the cosmic scale factor, from the past to the future, can be seen by 
following our curves from right to left. The circular dot on the vertical axis represents our 
"current time" (x = 0). We can clearly see that the only solution which shows a redshift in 
the past (values below the horizontal black dashed line, representing z = 0) is the red-solid 
curve, corresponding to k = —1. Therefore, the other two solutions, for k = +1, 0, are ruled 
out; only the k = — 1 solution will be considered in the following. 

Our preferred solution in Fig. [TJ (red-solid) also shows a blueshift region in the immediate 
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past of our current time, which in Sect. Ill II will be related directly to the Pioneer anomaly. 
This blueshift region is greatly exaggerated in the figure, since the different curves were 
plotted for 5 ~ 0.6, an unrealistically high value. We will show in the next sections that 5 
is positive and close to zero, resulting in a very small-sized blueshift region, compared to 
the overall size of the Universe. Similar plots can be obtained for the ratio R/Ro expressed 
in terms of the radial distance r (see 14| for details), which also suggest the existence of a 
blueshift region localized around the observer's position, i.e., the Earth could be surrounded 
by a natural blueshift region, extending at least over the solar system region. This might 
be the origin of the Pioneer anomaly. 2 



Before we proceed to analyze this possib 
more results obtained in our second paper 



e explanation for the anomaly, we recall a few 



151 ] . Since we have closed- form expressions for 



our scale factor R, in Eqs. (j4|) and (J5J), it is straightforward to obtain the Hubble parameter 

R 2 (t) 



[H(t) = R{t)/R{t)) and the deceleration parameter (q(t) = — azgop m — — fl C0 fl (*) ^ as a 



function of time or redshift z. For the k = — 1 case, we obtained [1 



H(t) = VW\c ( Sinh /"t C °t X l = ±VW\c Vil + f~[ 1 ^ 1 (7) 
\ cosh x — o smh x J (l + ^J 

/jx = ( coshx - £sinhx \ 2 _ 2 = (1 + zf ^ 

\ sinh x ~ 8 cosh x J (1 + z) 2 — (1 — 5 2 ) 

and, in particular for x ~ > or z — > 0: 



H(t ) = -±c ; H(z = 0) = ±jc (8) 
q(t ) = q(z = 0) = i - 2. 

The signs of the quantities in Eqs. ([7]) and (jHJ) can be explained by considering again the 
red-solid curve in Fig. [fl which represents the ratio R(x)/R{0), or equivalently R(t)/R(t ), 
over different cosmological epochs. This bell-shaped curve was plotted for a positive value 

2 Obviously, the Earth's observer is not located at any privileged position. The same cosmological evolution 
described by CC would be seen by any other observer in the Universe, provided that the local values of 
the cosmological parameters S and 7 are the same. In our previous work [3]) we have suggested 

that S might play the role of a universal time, so that for a certain value of this parameter the evolution 
of the Universe would look the same for any observer. In this way, conformal cosmology does not violate 
the cosmological principle, which postulates a homogeneous and isotropic Universe. 
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of 5 and shows a local blueshift area in the "past" evolution of the Universe, extending back 
to Xrs = arccosh [(1 + 5 2 ) /(l — 5 2 )] = 2 arctanh 5 (represented by the square point in Fig. [I]) 
or (t — t rs ) = J— arctanh 5, for the look-back time at which the redshift (rs) starts being 

W\k\c 



observed. The red curve has a maximum at Xmax = arctanh 5 or (to — t) mBX = i— arctanh 5 

W\k\c 



(we can also find (R(x) / R(0)) max = 1/x/l — S 2 or z min = y/1 — 5 2 — 1) and it is evidently 
symmetric around this point of maximum expansion of the Universe. 

Therefore, for each value of z, i.e., for each value of R(x)/R{0), we have two corresponding 



values of the Hubble parameter (except at the maximum, for z m ; n = a/1 — 5 2 — 1, where 
H = 0). The two related points on the curve, at the same redshift level, will have equal 
and opposite expansion rates. This yields the double sign in the previous expressions for H, 
when given as a function of z. This argument applies also to the z = case, corresponding to 
the current time to, at which H(t ) = —\c is negative, showing that the Universe is already 
in a contracting phase. 3 As discussed above, the same z = value can also refer to the time 
in the past (t rs ) at which we start observing the cosmological redshift, with H(t rs ) = +|c, a 
positive quantity. This analysis does not contradict the current astrophysical estimates of Hq 
as a positive quantity. They are based on redshift observations of light coming from galaxies 
at times in the past t < t rs , therefore what is denoted by H in standard cosmology should 
be actually indicated as H{t rs ) = +|c, again a positive quantity related to the expanding 
phase of the Universe. The same analysis can be done in terms of radial distances r. The 
blueshift region would extend from r = up to a distance given by: 

7 2 4 5 2 

r„ = -r/(l*l-^-) = - T 3a5. < 9 > 

where r rs is the distance at which we start observing the cosmological redshift. In general, 
the slope of the red-solid curve in Fig. [1] is related to the value of the Hubble parameter 
at that point, while its curvature is connected to the deceleration parameter, through the 
expressions given above. 

In particular, following Eq. flH]), the slope of the plot and its curvature at current time to 
are basically connected to our two fundamental parameters 7 and 5. In the next section we 



3 This is also a consequence of the signs of our conformal parameters, in particular the positive value of 7. 
Our estimate of 7 will be given in Sect. IIII| but we recall that Mannheim has independently evaluated 
7 as a small but positive quantity (jMann = 3.06 x 10 -30 cm" 1 ), by fitting rotational velocity curves for 



several spiral galaxies, using conformal gravity 13j. If 7 were to have a negative value, we would still be 
in an expanding phase of the Universe. 
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will show that the slope of the red-solid plot at to is closely related to the value of the Pioneer 
anomalous acceleration ap, which can therefore be used to determine 7. Similarly, the 
curvature of the plot at to wm be related to the rate of change of the anomalous acceleration 
(i.e., the "jerk" jp = ap) and will be used to determine the value of our other parameter 5. 

We conclude this section by noting that the values of our parameters (5 and 7) could 
be derived directly from standard cosmological observations, in view of Eq. (jSJ). Using the 
current best estimate of Hq = (72 ± 3) km s' 1 Mpc~ l 
we obtain: 



19| and the positive sign in Eq. (jSJ) 



7 



2H 



;i.56±0.06) x 10 



~ 28 cm' 1 . 



(10) 



The direct determination of 5 is more difficult, since the deceleration parameter q is not 
known explicitly. In Ref. 15j we based our analysis on recent luminosity data for type-la 
Supernovae, obtaining an estimate of 5 ~ 3.83 x 10~ 5 , but this analysis needs to be confirmed 
by further studies. 



III. THE PIONEER ANOMALY 



In the previous section we briefly reviewed our conformal cosmology and outlined the 
reasons why we consider the k = — 1 solution as a possible description of the evolution of 
the Universe. This solution can explain the observed cosmological redshift, but it requires 
the existence of a blueshift region in the immediate vicinity of our current spacetime position 
in the Universe. 

This could be a serious problem for our model, since we do not observe blueshift of 
nearby astrophysical objects except for the one caused by the peculiar velocities of nearby 
galaxies, presumably due to standard Doppler shift. However, as already mentioned in Sect. 



ocal region of blueshift might come from the analysis of the 

3, EL Q, Q, H, 0)- 



161, 20] 



HI experimental evidence of a 
Pioneer anomaly j^J, 

This is a small frequency drift (blueshift), observed analyzing the navigational data of 
the Pioneer 10-11 spacecraft, received from distances between 20 — 70 AU (astronomical 
units) from the Sun, while these spacecraft were exploring the outer solar system. This 
anomaly is usually reported as a positive rate of change of the signal frequency, Vp > 
(blueshift), resulting in a frequency drift of about 1.5 Hz every 8 years, or as an almost 
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constant sunward acceleration, ap < 0, or even as a "clock acceleration" a t = — < 0. More 
precisely (fl, £|): 



up = (5.99 ±0.01) x KT 9 s 
a P = -(8.74 ± 1.33) x 10~ 8 



(11) 



cm s 



a t 



ap 

c 



-(2.92 ±0.44) x 10~ 18 s' 1 



An attempt was made to detect such anomaly also in the radiometric data from other 
spacecraft traveling at the outskirts of the solar system, such as the Galileo and Ulysses mis- 
sions [20(. In the case of Galileo, the effects of solar radiation made such detection impossible, 
while for Ulysses a possible anomalous acceleration aui yS ses = — (12 ± 3) x 10 -8 cm/ s 2 was 
seen in the data. Other spacecraft, such as the New Horizons mission to Pluto, launched in 
2006, might provide new data in the near future. These discoveries prompted a complete 



re-analysis of a. 

n 

underway ([1], 



1 the historical navigational data of these space missions, which is curre ntly 
23 1, Q, 0, Q) and will be completed in the next few months 29]. 



This new analysis will try to determine additional characteristics of the anomaly, such as its 
precise direction, the possible temporal and spatial variations, its dependence on heliocentric 
or geocentric distance, etc. A future dedicated mission is also being proposed (3(3], 31], 



321 ]. [331 ]) to test directly this puzzling phenomenon. 

Currently, the origin and nature of this anomaly remains unexp 
of systematic errors have been considered ([lj], 0], 20], 21], 24], 



ained; all possible sources 



28], [34j ) but they cannot 



fully account for the observed effect. The current focus of conventional explanations of the 
anomaly seems to be the thermal recoil force, i.e., anisotropically emitted thermal radia- 
tion, originating from the spacecraft four radioisotope thermoelectric generators (RTGs), 
which can contribute significantly to the measured acceleration. The natural decay of the 
radioactive material in the RTGs, the aging of the thermocouples in the system and other 
effects, all contribute to the decrease of the total thermal power during the spacecraft life. 
This might explain the decrease over time of the measured Pioneer acceleration (in absolute 
value), i.e., the negative "jerk" ap < 0, already seen in the early Pioneer data ([lj], jl6 |. 

3). 

Although the anomaly can be caused by standard physical effects, we will try in the 
following to explain its origin by using the cosmological model outlined in the previous 
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section. The phenomenology of the Pioneer anomaly is related to a complex exchange of 
radiometric signals between the tracking stations on Earth (of the Deep Space Network - 
DSN) and the spacecraft, using S-band Doppler frequencies (1.55 — 5.20 GHz). Typically, 
an uplink signal is sent from the DSN to the spacecraft at a frequency of 2.11 GHz, based on 
a very stable hydrogen maser system, then an S-band transponder onboard the spacecraft 
applies an exact and fixed turn-around ratio of 240/221 to the uplink signal, so that the 
Pioneer returns a downlink signal at a slightly different frequency of about 2.29 GHz, to 
avoid interference with the uplink one. 

This procedure is known as a two-way Doppler coherent mode and allows for very pre- 
cise tracking of the spacecraft, since the returning signal is directly compared to the original 
one. On the contrary, a one-way Doppler signal (with a fixed signal source on the spacecraft, 
whose frequency cannot be monitored for accuracy) is less effective. This type of tracking 
system added to the propulsion and navigational characteristics of the Pioneer spaceship 
(especially the presence of a spin-stabilization system) resulted in a very good accelera- 
tion sensitivity of about 10~ 8 cm/s 2 , once the influence of solar radiation pressure can be 
neglected (for distances > 20 AU from the Sun). 

The DSN station acquires the downlink signal after a time delay ranging from a few 
minutes to some hours, depending on the distance involved, and compares it to the reference 
frequency to determine the Doppler shift due to the actual motion of the spacecraft. The 
navigational software can also model with great precision the expected frequency of the 
signal returned from the Pioneer, which should coincide with the one observed on Earth. 
As already mentioned, a discrepancy was found, corresponding to the values in Eq. (TTTj) . 
whose origin cannot be traced to any systematic effect due to either the performance of the 
spacecraft or the theoretical modeling of its navigation. 

The Pioneer anomaly was first reported [3], j4|) as an almost constant value of the 
anomalous acceleration, with temporal and space variation of ap within 10%, over a range of 
heliocentric distances ~ 20 — 70 AU, and possibly at even closer distances < 10 AU, so that 
we will concentrate first on the average value of ap and later on its variation with time and 
distance. In our view, the Pioneer phenomenology represents the most basic experiment we 
could perform in order to check if the cosmic evolution is really affecting the frequency of 
electromagnetic radiation emitted and observed at different spacetime locations, following 
Eqs. © and ©. 
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In the standard analysis of the Pioneer anomaly, the signal coming back to Earth is 
affected by the relativistic Doppler effect. Following this model, z/ mo d will be the frequency 
of the expected signal and will be related to the signal reference frequency v re f = 2.11 GHz 
(for the uplink signal in a two-way system) by the standard relativistic Doppler formula (see 
Eq. 2.2.2 in |35|): 



^ = ^^-1-7, (12) 

where v r is the spacecraft radial velocity and the approximation on the right-hand side holds 
to first order in v r /c. 

Since we have a two-way system, the Doppler shift involved is actually double, so we can 
use the previous equation but with v r = 2v mo d(t'), where v mo d(t') is the expected velocity of 
the spacecraft, according to the theoretical navigation model, at time t', when the spaceship 
receives and immediately re-transmits the signal. We use here a time variable t' which can 
be simply considered the elapsed time since the spacecraft launch (f — at r — 0) and 
then later we will simply identify t' with our cosmological look-back time (to —t) in Eq. ([61). 
With this radial velocity, Eq. ( [121) to first order in v r /c becomes: 



^mod(^ ) 



'ref 



_ 2v mod (t') 



(13) 



c 

and this frequency is expected to be observed with high precision, due to the reported 
excellent navigational control of the spacecraft. 

On the contrary, a different frequency is observed, v bs(t') > v m od(t'), involving an addi- 
tional unexplained blueshift: this is the Pioneer anomaly. Following Eq. (jlip . the frequency 
difference is reported as: 



Ai/(0 = v obs (t') - zVod(0 - 2tV P 
up = 5.99 x 1(T 9 s~ 2 (one-way) 



(14) 



where the factor of two in the first line of the previous equation is due to the two-way system. 
We also remark here that several of the cited references adopt a rather confusing "DSN sign 
convention" for the frequency difference in Eq. (j!4p (see jj, 4j, 23| and Ref. (38) of |2j), 
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resulting in a change of sign in most of their equations. We prefer to use here our definition 
of Av as given in the previous equation. 

The anomalous acceleration ap is introduced as an alternative way of describing the 
effect, although in our view it does not correspond to a real spacecraft acceleration. As in 
Eq. (II 3p . we can write the observed frequency to first order in v r /c as: 

;1 _ 2v obs (t') 



V obs (t') ^ Vref , (15) 

where the "observed" velocity of the spacecraft refers to the time of interest t' . Combining 
together the last three equations we can write the frequency difference as: 

Az,(f) = -2^L [Vobs{t >) _ Vmod{t ')] = -2^lAv(t'). (16) 
c c 

These frequency differences Av (also called frequency residuals in the literature cited) are 
therefore equivalent to the corresponding velocity residuals (Av = v Q b s — f mod) and they are 
usually plotted as a function of the elapsed time t', showing an almost linear increase with 
time of these residuals, which is the essence of the Pioneer anomaly (see for example Fig. 
5.2 in Ref. [l|). The Pioneer anomalous acceleration can be defined as the rate of change of 
the velocity residuals, related to the corresponding rate of change of the frequency residuals, 
in view of Eq. §Uo§. Therefore, if we define a P = ~ [Av(f + At') - Av(t')\ /At', the 
Pioneer acceleration can be related to the frequency differences 

c d(Au) c Au(t' + At') - Av(t') 

a P = ^ — - ~ - 1 — , 17) 

2v ref dt' 2u ref At' v ; 

which are more significant quantities in our analysis. We will assume that these frequency 

differences are intrinsically due to the different locations of the spacecraft (at position r) 

and of the Earth's observer (at r = 0). Therefore, we identify the reference frequency v re j in 

Eq. (fT3]) with v(0) and the similar quantity u re f in Eq. ffl5l) with z/(r). Then, we subtract 

Eq. (QJ) from Eq. ([151): 



Av(t') = 2 [v(r) 



v 



1 - 



2 m 



ref 



~R{r) 
W) 



2v(t') 

c 

2v re f 



2 \v{r) - z/(0)] = 2i/(0) 



u(r) 

MO) 



(18) 



m 

R(0) 



2v, 



ref 
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where the common factor of two in all the parts of the previous equation was added again 
because of the two-way effect, which has to be included also in our gravitational blueshift 
model. The velocities v mo d(t') and v b s (t') from Eqs. (ITS]) and (TO)]) are assumed to be the 



same, so that the common factor 



1 



2v(t') 



1 is close to unity and can be neglected, 



since the average Pioneer speed is vp — 12.8 km/s <C c |36[]. We also identified z/(0) with 
the Earth reference frequency v re f and used our fundamental Eq. (jlj) and Eq. k < 
case, to first order in x — \fW\ c ^' Similarly, we have: Au(t' + At') ~ 2u re f [| (t' + At')] 
so that Eq. (fT7|) simplifies as follows: 

a P = a P (t ) = -|c 2 = cH(t ), (19) 

in view also of our evaluation of H(to) = — (a negative quantity) from Eq. (JE]). 

This result immediately explains the often cited "numerical coincidence," i.e., the simple 
relation \ap\ ~ cH Q between the Pioneer acceleration and the standard (positive) Hubble 
constant, with the correct negative sign for both quantities in Eq. (|T9|) . in view of our 
previous discussion of the sign of H(t ) < 0. Eq. f fT9"j) can also be used to determine 7 and 
H (as a positive quantity), using the reported value of ap from Eq. ( TTTi) : 

7 = 7 (t ) = -%a P = (1-94 ± 0.30) x 10~ 28 cm- 1 , (20) 
cr 

H = (90.0 ± 13.7) km s' 1 Mpc' 1 . 

The value of 7 (considered measured at the current time to, even if the Pioneer data are 
a few years old) is close to our first direct estimate in Eq. (11 01) and the corresponding value 
of the Hubble constant is close to the value of standard cosmology. We remark here again 
that our model fully explains the reason of this "numerical coincidence" and provides also 
the correct signs for all the quantities involved. 5 



4 The elapsed time t' for the Pioneer spacecraft missions is of the order of a few years (lyr — 3.156 x 10 7 s); 
we can assume yj\k\ ~ 7 ~ 1CT 28 - 1CT 30 cmT 1 , therefore ^/\k\ct' — 1CT 10 - 1CT 12 < 1. 

5 The numerical "coincidence" between the Hubble constant and the value of the Pioneer acceleration 
ap divided by c, was noticed immediately after the discovery of the Pioneer effect and prompted many 
speculations and different explanations. This coincidence is even more striking if one uses the value cited 
in Ref. Q as the experimental value for Pioneer 10 data before systematics, ap = —7.84 x 10~ 8 cm s -2 , 
thus obtaining H — 80.7 km s" 1 Mpc^ 1 and 70 = 1.74 x 10 -28 cm -1 . 



14 



Following Eqs. (|T7|) -( TT9|) and the related discussion we can generalize our expression of 
the Pioneer acceleration, SIS db function of time t'\ 



d[R(t')/R(0)} _ 27 [sinhx-^coshx] (oU 
dt 2o [cosh x — 8 sinh x) 

with x — \/\k\ct' = ^ct'. 6 In particular, by using the previous equation and taking another 

time derivative, it is easy to derive the "jerk" j P = dap/dt' and its value in the limit for 

f -> 0: 

JP=Jp(to) = c 3 (^) 2 (l-25 2 ), (22) 

expressed in terms of our fundamental parameters 7 and 5. The current value of jp in the 
last equation is positive (for small values of S), but the Pioneer acceleration, as in Eq. (TlTI) or 
Eq. ffT9]) . is considered negative in this paper so that a positive jerk means that the absolute 
value of dp will decrease for increasing times or radial distances, which is indeed shown in 
the early Pioneer data, as it was already mentioned at the beginning of this section. 



In Fig. [2] we illustrate the early Pioneer 10/11 data, as originally reported in Ref. 16], 
where the absolute value of the Pioneer acceleration \ap\ is plotted as a function of the 
radial heliocentric distance in AU. The red-dashed horizontal line and the green-dotted 
lines represent respectively the value of \ap\ and the related 1-sigma error range quoted in 
Eq. (fTT|) . The first three data points for Pioneer 11, at smaller distances, lie outside the 
considered error range probably because the anomalous acceleration was masked by solar 
radiation or other effects. We will not include these first three data points in our subsequent 
discussion. We will concentrate our analysis on either just the Pioneer 10 data points, or the 
combination of data points for both spacecraft, but within the 1-sigma error range ("error 
range data" in the following). 

These two sets of data clearly show a possible decrease of the Pioneer anomaly (in the 
absolute value \ap\) with increasing heliocentric distance. The black (dash-dotted) line and 
the blue- solid line in the same figure represent linear fits for the Pioneer 10 and the error 
range data respectively, both of them indicating a decrease of \ap\. 



6 Although t' is the elapsed time since the spacecraft launch, it is treated here as equivalent to a look-back 
time (to — t) because the Pioneer is moving toward increasing distances r, therefore corresponding to 
increased look-back times in our original redshift interpretation. 
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FIG. 2: Early data for Pioneer 10/11 acceleration as a function of heliocentric distance. The 
average value of the anomalous acceleration is indicated in red-dashed, together with its error 
range (green-dotted). We also show linear fits of the data, which allow for the determination of 
our cosmological parameters 7 and 5. 

If our conformal cosmology is the origin of the Pioneer anomaly, and not the thermal 
recoil force mentioned at the beginning of this section, our "jerk" equation (122]) will explain 
the decrease of \ap\ and can also be used to determine our second parameter 5. 

We computed the slopes of our two linear fits in Fig. [2] and used them as (positive) values 
of jp in Eq. (I2"2"j) . together with the 7 value from Eq. (I20p . 7 Solving Eq. ( 12 2 p for 5, we 

7 The radial distances r of the data plotted in Fig. [5] were converted into elapsed times t', by using a 
simple approximation: r ~ vpt', where vp is the average Pioneer speed. From the original data (available 
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obtain: 



j P = (3.85 ± 1.88) x 1(T 17 cm s~ 3 (Pioneer 10 data) (23) 

5 = (8.12 ±2.35) x 1(T 5 
jp = (1.37 ± 0.95) x 10~ 17 cm s~ 3 (error range data) 

5 = (1.36 ±0.52) x 10~ 4 



15| (S = 3.83 x 10" 5 ), 



and these values for S are very close to the one we obtained in Ref. 
which was based solely on the analysis of type-la Supernovae data. 

Another type of analysis is illustrated in Fig. [3J The Pioneer 10/11 data, the standard 
value of \ap\ and the related error range are the same as in the previous figure, but this 
time we used the generalized expression of ap in Eq. ( l2Tj) to fit the data within the error 
range. We allowed both quantities 7 and 5 to be free parameters in our fitting procedure 
and we converted the elapsed time t' in Eq. ( 121]) into the radial distance r by using the 
approximation r ~ vpt', where vp is the average Pioneer speed, as it was done also for the 
data in the previous figure. The radial distance r should be more properly identified with 
the geocentric distance of the spacecraft, rather than the heliocentric one, since r should be 
the distance from the Earth observer. We also performed fits using the geocentric distance, 
but the results were very similar to those obtained by using heliocentric distances, so we 
will not include them in the following analysis. 

Again, in Fig. [3] we used the expression in Eq. ( 12 ip to fit the data, although the fitting 
curves appear almost as straight lines in this figure. The first conformal cosmology fit, 
illustrated by the black (dash-dotted) curve, was obtained by using only the Pioneer 10 data 
and yielded the following values of the parameters: 



5 = (9.19 ± 1.53) x 10~ 5 (Pioneer 10 data) (24) 
7 = (2.20 ± 0.18) x 10~ 28 cm- 1 . 

The second fit (blue-solid curve) was obtained by using all the data within the error range 
(again omitting the first three Pioneer 11 data points) and produced the following results: 



from the NASA website at: http://cohoweb.gsfc.nasa.gov/helios/) we estimated: vp W ~ 12.96 km/s 



vpn ~ 11.42 km/s, and used an average vp ~ 12.19 km/s when combining data for both spacecraft. 
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FIG. 3: Early data for Pioneer 10/11 acceleration as a function of heliocentric distance. The 
average value of the anomalous acceleration is indicated in red-dashed, together with its error 
range (green-dotted). We also show full conformal cosmology fits of the data, which allow for a 
better determination of our cosmological parameters 7 and 5. 



5 = (1.38 ± 0.43) x 10" 4 (error range data) (25) 
7 = (1.97 ± 0.08) x 10~ 28 cm- 1 . 

Comparing the results in the last two equations with those for S in Eq. f[2"3"j) . obtained 
with a fixed 7 = 1.94 x 10~ 28 cm -1 as in Eq. ( 1201) . we can see that all the values of our 
parameters are in agreement. In particular, from the different analyses we consistently 
obtain 7 ~ 1.9 — 2.2 x 10~ 28 cm -1 and 5 ~ 10~ 4 — 10 -5 , where the different values depend 
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on the Pioneer data being used. As already remarked, the values for 5 quoted above are 



also very close to the one we obtained in Ref. [15] (So = 3.83 x 10~ 5 ), based on type-la 
Supernovae data. In the next section we will discuss our results and compare them to the 
current limits of standard gravity in the solar system. 

IV. DISCUSSION OF OUR RESULTS AND CONCLUSIONS 

In the previous sections we discussed how conformal cosmology provides a natural expla- 
nation for the Pioneer anomalous acceleration, in both magnitude and direction (i.e., the 
negative sign of the radial acceleration). We also explained the "numerical coincidence," con- 
necting ap with the Hubble constant, and the observed decrease with heliocentric distance 
of \ap\, related to the Pioneer jerk jp. Although the Pioneer data are still not very accurate, 
our analysis consistently indicated that our conformal parameters are approximately given 
by 7 ~ 1(T 28 cmT 1 and 5 ~ 1(T 4 - 1(T 5 (see Eqs. (PHJ» , (J2DJ and (P])-(j25])). In this final 
section we will discuss the implications of the values of our parameters in relation to other 
studies in the field. 

We first remark that a new analysis of rotational velocity data for spiral galaxies, based 
on conformal gravity, has recently appeared js3], improving the original work on the subject 



171 ]. 18j). This new study uses the full line element of conformal gravity in Eqs. (|T])-( 



including the effects of the quadratic term —nr 2 , which were previously neglected, thus 
obtaining a global gravitational potential V g i b a i(r) = ^c 2 r — ^c 2 r 2 of cosmological origin. In 
addition to this, a local gravitational potential V/ oca /(r) is obtained by integrating over the 
visible galactic mass distribution a gravitational potential per unit solar mass of the form 
V*(r) = —G^f 1 + ^c 2 r. The two potentials, global and local, are then combined together 
to model the rotational motion of galaxies. The fits to galactic rotation data [37J, performed 
without any dark matter contribution, show a remarkable success of conformal gravity, even 
at the largest distances from the galactic centers, where the quadratic term — nr 2 becomes 
important and comparable to the linear term 77*. Mannheim and collaborators (3^ were then 
able to determine the values of the global universal parameters as jMann — 3.06 x 10~ 30 cmT 1 
and KMann = 9.54 x 10~ 54 cm~ 2 . The related terms of the global gravitational potential were 
associated respectively to the cosmological background and to cosmological inhomogeneities. 
The local parameter 7* was also evaluated as 7* = 5.42 x 10~ 41 cm -1 . 
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The values of the dimensionful parameters 7 and k obtained through this analysis of 
galactic rotation curves are somewhat different from our values, reported in this paper or in 
our previous work [jjj] (70 = 1.94 x 10 -28 cm -1 and k = 6.42 x 10~ 48 cm -2 ). This difference 
could be due, as we explained in Ref. jisj, to a possible redefinition of the luminosity 
distance and other distance indicators, which might affect even the radial distances (from 
the galactic centers) which are employed in the galactic rotation analysis. 

However, it is instructive to compute the dimensionless S parameter, using Mannheim's 
values 7 Mann and i^Mann, because this dimensionless constant should not be affected by 
a revision of the cosmo logical distances. As explained at the beginning of Sect. [TTJ the 
parameters k, 7, k and 5 are related through k = — \ — k and also 5 = (7/2 y/\k\) so that 
we obtain: 



k M ann ^ ~^Mann = "9.54 X 10 54 CTU 2 (26) 
S M ann = 4.95 X 10~ 4 . 

Therefore, the conformal gravity analysis by Mannheim and collaborators suggests a k = 
||t = — 1 Universe, consistent with our cosmological model and also a value of SMann = 
4.95 x 10~ 4 , close to our quoted values of 5 ~ 10~ 4 — 10~ 5 . 

Conformal gravity considers local gravitational effects as being due to the local potential 
Viocai, or simply to the potential V*(r) = —G^- + \(?r for our solar system. Since the 
value of the local constant 7* ~ 10~ 41 cm -1 is very small, the modifications to standard 
dynamics of the solar system are negligible. 8 Therefore, conformal gravity is not in any way 
in contradiction with the very stringent limits on alternative gravity theories imposed by 



studies of planetary ephemerides, or other solar system observations (|5|, [38|], [39|, [40| ) . 

As for our analysis of the Pioneer anomaly, we used the reported values of the anomalous 
acceleration ap to determine the cosmological parameters, simply because such was the 
way these data were reported in the literature cited. However, it should be clear from the 
discussion in Sect. IIHI that we explain the Pioneer anomaly in terms of our cosmological- 



For example, the ratio between the conformal gravitational potential -^-c 2 r and the standard Newtonian 
term G^p- at a heliocentric distance of 1 AU is ~ 10~ 20 , while at a distance of 100 AU (outer solar 
system) the same ratio is ~ 10 -16 . Therefore, the "conformal gravity force" is negligible, compared to 
the standard Newtonian one, over the whole solar system region. 
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gravitational blueshift, based on the global values of the parameters 7, k and 5. In this view, 
there is no real dynamic acceleration of the Pioneer spacecraft (or of any other object in the 
solar system) oriented toward the Sun, due to some new gravitational force or modification 
of existing gravity, except for the tiny corrections coming from local conformal gravity 
mentioned above. In fact, in our analysis we assume that there is no difference between the 
two velocities v mo d(t') and v b s (t') in Eqs. ( !T3|) and (|T5|) . therefore the anomalous acceleration 
defined as ap = m — [Av(t' + At') — Av(t')] /At' with Av = v bc— v mo d is actually zero. 

In this way we also overcome the objection, reported in Ref. |4|, that "the anomalous 
acceleration is too large to have gone undetected in planetary orbits, particularly for Earth 
and Mars," since "NASA's Viking mission provided radio-ranging measurements 41] to an 
accuracy of about 12 m," which should have shown the effect of the anomalous acceleration 
on the orbits of these two planets. 

In our view, precision ranging measurements with radio signals or lasers, based on the 
round-trip travel time from Earth to other bodies in the solar system, would not show any 
anomalous effect because the speed of light is not affected by our cosmological model and 
the corrections to the dynamics of the solar system due to conformal gravity are negligible. 

On the contrary, we would observe an effect similar to the anomalous acceleration for 
a spacecraft, a planet, or any other object in the solar system, if we were to study its 
motion through Doppler frequency ranging, because of the intrinsic differences in frequency 
or wavelength for light emitted at different spacetime positions, due to our cosmological 
model. 

The size of the local blueshift region, which in our model is responsible for the frequency 
differences, can be easily estimated by using Eq. ([9]) and the values of our parameters. 
For example, using the values from our conformal cosmology fits in Eqs. (124]) and fl25|) . 
we obtain r rs ~ 50 — 126 pc, corresponding to a distance comparable to the one between 
Earth and the nearest bright stars (which is about 15 — 30 pc). This blueshift region would 
extend well beyond the solar system, but would cover a small portion of our galaxy, since 

TMilkyWay — 14.6 kpC. 

The maximum blueshift effect would be seen at r = \r rs ~ 25 — 63 pc and would 



correspond to a z min = y/1 — 5 2 — 1 ~ — 10~ 8 , a very small value. Therefore, the blueshift 
region and the related effects are so small that they cannot be practically observed in the 
radiation spectrum of stars or other radiation emitting objects within this region. These 
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effects are only small corrections to the Doppler signals coming from the Pioneer or other 
similar spacecraft. 

Finally, we want to compare our estimates of the rate of change of the anomalous accel- 
eration (i.e., the jerk jp) with those presented by independent verifications of the Pioneer 
anomaly (see review in Ref. The first of these studies was performed by Markwardt 



361 ]. who reviewed data for Pioneer 10 and reported apio = —(7.70 ± 0.02) x 10~ 8 cm/s 2 , 
with jpiQ < 0.18 x 10~ 8 cm/s 2 /year = 5.70 x 10" 17 cm/s 3 . 9 Using Markwardt values in 
Eqs. (TT9|) and (1221) we obtain ^Mark = 1-71 X 10~ 28 cm" 1 and duark = 5.89 x 10~ 5 , consistent 
with our values in Eq. (I20p and Eqs. (j2"3"l) - (j2"4"j) for Pioneer 10. 



The second independent study was done by Toth 42| and reported results separately for 
the two spacecraft. From Toth's results for Pioneer 10 (apio = — (10.96±0.89) x 10" 8 cm/s 2 , 
jpio = (0.21 ±0.04) x 10 -8 cm/s 2 /year = 6.65 x 10~ 17 cm/s 3 ) we compute jTothPio = 2.44 x 
10~ 28 cm -1 and SrothPio = 7.76 x 10~ 5 . Using instead Toth's results for Pioneer 11 (apn = 
-(9.40 ±1.12) x KT 8 cm/s 2 , j Pn = (0.34 ±0.12) x 10~ 8 cm/s 2 /year = 1.08 x 10" 16 cm/s 3 ) 
we obtain ^TothPU = 2.09 x 10~ 28 cm" 1 and SrothPii = 5.23 x 10~ 5 , and all these results are 
also consistent with those discussed in Sect. IIHI 

In conclusion, the detailed analysis of the Pioneer anomaly presented in this work has 
indicated that our conformal cosmology might be the origin of this effect, while conformal 
gravity alone cannot account for the anomalous acceleration of the spacecraft. If our analysis 
is correct, it explains naturally the numerical coincidence between the Pioneer acceleration 
and the Hubble constant, including the signs of these quantities. In addition, we confirm our 
previous evaluations of the cosmological parameters, 70 = (1.94 ± 0.30) x 10~ 28 cm -1 and 
S = 3.83 x 10" 5 , also in agreement with independent evaluations. Further studies will be 
needed when the re-analysis of all the historical navigational data for the Pioneer spacecraft 
will be completed by S. Turyshev and collaborators and new data will be publicly available. 



9 We prefer to report here, as also done in the rest of the paper, the anomalous acceleration spasa negative 
quantity and the related jerk jp as a positive quantity. Some of the papers in the literature adopt the 
opposite sign convention, which might generate some confusion. 
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